0QIV-C1V  BOO  TUf,3U| 


U.S.  ARMY  TANK-AUTOMOTIVE  COMMAND 
RESEARCH  AND  DEVELOPMENT  CENTER 
Warren,  Michigan  48090 


NOTICES 


This  report  is  not  to  be  construed  as  an  official  Department  of  the  Army 
position. 


Mention  of  any  trade  names  or  manufacturers  in  this  report  shall  not  be 
construed  as  an  official  indorsement  or  approval  of  such  products  or 
companies  by  the  US  Government. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the 
originator. 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECjJBl.TY  CLASSIFICATION 

Unclassified 

lb., RESTRICTIVE  MARKINGS 

hone 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION /.AVAILABILITY  OF  REPORT  , .  ... 

Approved  for  public  release  distribution 

2b.  DECLASSIFICATION /IX/WNGRAOING  SCHEDULE 

unlimited 

4.  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 

5.  MONITORING  ORGANIZATION  REPORT  NUM8ER($) 

TACOM  T.R.  #13004 

8a.  NAME  OF  PERFORMING  ORGANIZATION 

Michigan  State  University 

6b  OFFICE  SYMBOL 
(If  applicable) 

7a.  NAME  OF  MONITORING  ORGANIZATION 

TACOM  AMSTA-RCKM 

8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 


7b.  ADDRESS  (Oty,  State.  and  ZIP  Code} 

TACOM 

Attn:  AMSTA-RCKM 
Warren,  MI  48090 


8b.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Of  applicable) 


Same  as  7 


8c  ADCRESS  (Oty,  SUM.  and  ZIP  Cot*) 


Contract  DAAE07-82-K-4Q02 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


1 1  TITLE  ( Induct *  SKunty  Otsufkttfon) 


Theoretical  and  Experimental  Investigation  of  Mechanically  Fastened  Composites 


"  2  PERSONAL  AUTHQR(S) 

Cloud,  Gary;  Sikarskie,  David;  Mahajerin,  Enayat;  and  Herrera,  Pedro 


14.  DAT^OF^REPORT  (Year,  Month.  Day)  115.  PAGE  COUNT 


13a.  TYPE  OF  REPORT 

Final 


COSATt  CODES 


SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  idantify  by  block  number) 


19.  ABSTRACT  (Continue  on  reverie  if  necessary  and  idantify  by  block  number) 

The  mechanics  of  fasteners  in  composites  were  studied  1 r  a  combined  experimental - 
theoretical  research  program.  The  objectives  were  to  gain  fundamental  Insight  into  the 
stress-strain  field  near  pin-type  fasteners  and  to  provide  guidance  to  designers  responsible 
for  the  selection  and  sizing  of  fasteners.  The  primary  experimental  method  utilizes  Moire 
Interference  with  optical  Fourier  processing  of  grid  photoplates.  A  new  technique  using 
optical  Interference  to  generate  gratings  In  three  directions  was  developed  to  gain  a  factor 
of  10  In  sensitivity.  The  constitutive  properties  of  the  material  studied  were  measured. 

For  the  analytical  work,  a  boundary  element  method  (BEM)  was  developed,  a  compact  and 
efficient  computer  code  written,  and  the  method  compared  with  the  finite-element  method, 

~  (F£M)»-  cor  the  problems  Investigated,  the  BEM  Is  more  efficient.  The  material  used  was 
fiber  g»ass-epoxy  laminate  with  woven  fibers.  The  analytical  and  experimental  forces  were 
brought  to  bear  on  the  problem  of  a  loaded  pin  snugly  fit  In  a  hole.  Results  from  the  two 
approaches  agreed  well  for  the  specific  composite. 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

,  CSuNCLASSIFIEOAJNLIMITEO  □  SAME  AS  RFT.  □  OTIC  USERS 


22a.  I  AME  OF  RESPONSIBLE  INDIVIDUAL 

Avery  H.  Fisher 


00  FORM  1473, 84  MAR  83  APR  Wition  may  b«  used  until  •xhaustad.  security  classification  OF  this  page 

AH  othar  aditiom  are  obwlttt.  '  —  —  " 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


Include  Ana  Code)  22c.  OFFICE  SYMBOL 

9  AMSTA-RSC 


rtrr 


c  j  A)  [(  tifrf i  . i- 


PREFACE 


The  research  described  in  this  report  was  conducted  at  Michigan  State 
University  In  the  Department  of  Metallurgy,  Mechanics  and  Materials 
Science  and  in  the  Center  for  Composite  Materials  and  Structures. 
Investigators  were  Dr.  Gary  Cloud,  Professor;  Dr.  David  Silcarskle, 
Chairman;  Dr.  Enayat  Mahajerin,  Post-doctoral  Student;  and  Mr.  Pedro 
Herrera,  Doctoral  Candidate.  Manuscript  preparation  was  by  Ms.  Arlene 
Klingbiel. 


DT1C 

ELECT E 


MAR  6  1985 


B 


§ 

-X 


I 


This  page  left  blank  Intentionally. 


I 


j •  *  ■  *  i 

*,.•/  •/  v 


V  V  V  >  1  •  •  -  •••  V  ^  i 


•>  *>  ■  «vv  ‘ 

A  .*•  >  /•'.*« 


#  f* 


TABLE  OF  CONTENTS 


Section  Page 

1.0.  INTRODUCTION  .  .  .  ..........  11 

1.1.  Background  . 11 

1.2.  Project  Summary . .  11 

2.0.  ANALYTICAL  DEVELOPMENT,  RESULTS  AND  COMPARISON  ....  .  12 

2.1.  Problem  Statement  .  12 

2.2.  8EM  Formulation . . . .  :  13 

2.3.  Example  Problems . .  .  ......  .  .  18 

2.4.  Comparison  of  BEM  and  FEM  ....  .  ........  18 

3.0.  EXPERIMENTAL  METHODS  AND  RESULTS . ...... . 27 

3.1  Introduction  .  .  . . 27 

3.2.  Specimen  Fabrication  and  Material  Specification  .  .  . . 28 

3.3.  Specimen  Gratings  .  .  ..........  28 

3.4.  Printing  Grating  Onto  Specimen  . .  .  .  31 

3.5.  Grid  Photography  . 31 

3.6.  Summary  of  Optical  Processing  to  Create  Fringe  Patterns  ....  34 

3.7.  Experimental  Results  . 35 

3.8.  Experimental  Evaluation  of  Elastic  Constants  for  Composite  .  .  46 

3.9.  Experimental  Results  and  Discussion  .  .  .  .47 

3.10.  High-sensitivity  Interferometric  Moire  Technique  .  49 

4.0.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  NUMERICAL  RESULTS  .....  52 

REFERENCES  .  ............  54 

APPENDIX  A  FUNDAMENTAL  SOLUTIONS  .  .  . .  ........  A-l 

APPENDIX  B  THE  COMPUTER  PROGRAM  "BEM"  .  . .  B-l 

APPENDIX  C  SIMPLE  OPTICAL  PROCESSING  OF  MOIRE-GRATING  PHOTOGRAPHS  .  C-l 
DISTRIBUTION  LIST . .  .  .  .......  D-l 


5 


This  page  left  blank  Intentionally 


LIST  OF  FIGURES 


Figure  Page 

2-1.  Geometry  for  the  Boundary  Element  Method . 14 

2-2.  Geometry  for  Computation  of  a*8Y's  16 

2-3.  Double-lap  Mechanical  Joint . 19 

2-4.  Geometry  for  the  Example  Problem .  20 

2-5.  Schematic  of  Boundary  Element  Subdivisions. .  .  .  .  21 

2- 6.  FEM  and  BEM  Subdivision  Schemes  for  the  Mechanical  Jcint 

Problem . . .  . . _ ...  26 

3- 1.  Dimensions  of  the  Test  Coupon . 29 

3-2.  Grating  Production  Process  . 32 

3-3.  Set-up  for  Grating  Photography  .  33 

3-4.  Photograph  of  Moire  Pattern  Showing  Displacements  Perpendicular 

to  Load  Line  for  Zero  Load  with  Pitch  Mismatch .  36 

3-5.  Photograph  of  Moire  Pattern  Showing  Displacements  Perpendicular 

to  Load  Line  for  400  lbs.  Load  with  Pitch  Mismatch .  37 

3-6.  Photograph  of  Moire  Pattern  Showing  Displacements  Parallel 

to  Load  Line  for  Zero  Load  with  Pitch  Mismatch  .  38 


3-7.  Photograph  of  Moire  Pattern  Showing  Displacements' Parallel 

to  Load  Line  for  400  lb.  Load  with  Pitch  Mismatch  .  .  i  .  .  .  39 


3-8.  Strain  Perpendicular  to  Load  Axis  (ex)  Along  Several  Lines 

In  the  Ligaments  of  Fastener  Hole . .  40 

3-9.  Strain  Perpendicular  to  Load  Axis  (ex)  Along  Several  Lines 

In  Bearing  Region  Near  Hole . .  41 

3-10.  Strain  Parallel  to  Load  Axis  (ey)  Along  Several  Lines 

Encompassing  a  Ligament  Region  .  .  42 

3-11.  Strain  Parallel  to  Load  Axis  (ey)  Along  Several  Lines  in 

and  Opposite  to  the  Bearing  Region  ....■■ . .  43 

3-12.  Specimen  Showing  Location  of  Strain  Gages  ..........  44 


3-13,.  Typical  Tensile  Stress-strain  Plots  for  Determination  of 
Mechanical  Properties  of  Orthutropic  Composite  used  in  this 
Project . . . 

3-14.  Distribution  of  Tensile  and  Compressive  Strain  (ey)  Near 

the  Fastener  Hole . . . 

3-15.  High  Sensitivity  Moire  Fringes  Showing  Displacements 

Perpendicular  to  the  Direction  of  Loading  for  Mo  Load  with 
Rotational  Mismatch . . . . 

3- 16.  High  Sensitivity  Moire  Fringes  Showing  Displacements 

Perpendicular  to  the  Direction  of  Loading  for  a  130  lb.  Load 
with  Rotational  Mismatch  . 

4- 1.  Comparison  of  Stress  Concentration  Factors  Obtained  by 

Numerical  and  Experimental  Techniques  .  .  .  . . .  . 


LIST  OF  TABLES 


Table  Page 

2-1.  Stresses  along  the  llgamer.t  'Mne  ab  for  an  orthotropic 
composite.  Uniaxial  Tension,  d/w  *  .5,  M  *  60  (boundary 
subdivisions)  d  3  0.5,  e  *  C.5,  h  *  1.5,  w  *  1.0 . .  22 

2-2.  Stresses  along  the  ligament  line  ab  for  an  orthotropic 

composite.  Cosine  distribution  of  tractions  on  the  upper 

half  of  the  hole,  d/w  *  .5,  N  *  60  (boundary  subdivisions) 

d  *  0.5,  e  »  0.5,  h  *  1.5,  w  *  1.0 . .  23 

2-3.  Stresses  along  the  ligament  line  ab  for  an  orthotropic 

composite.  Cosine  distribution  of  tractions  on  the  upper 

half  of  the  hole,  d/w  '•*  .2,  N  »  60  (boundary  subdivisions) 

d  «  0.2,  e  »  0.5.  h  *  1.5,  w  -  1.0 . . . .  24 

2-4.  Stresses  along  the  ligament  line  ab  for  an  orthotropic 
composite.  Fixed  pin  case,  d/w  *  .5,  N  *  60  (boundary 
subdivisions)  d  *  0.5,  e  ■  0.5,  h  *  1.5,  w  *  1.0 .  25 


9 


v- 


* 

3 

■s- 


'  v’ 


£r) 

"ja 


1.0.  INTRODUCTION 
1.1.  Background 


There  has  been  a  dramatic  Increase  1r.  che  use  of  composite  materials.  The 
center  of  gravity  of  this  growth  has  been  In  aerospace  and  defense-related 
Industries.  As  composite  material  developments  continue  to  Increase, 
resulting  In  decreases  In  cost,  use  of  these  materials  will  naturally 
spread  to  other  Industries. 

One  of  the  attractive  features  of  composites  as  opposed  to  metals  Is  their 
ability  to  be  molded  Into  complex  "final  shape."  Important  advantages  of 
this  approach  include  requiring  a  minimum  number  of  fasteners  as  well  as 
being  able  to  locate  the  necessary  fasteners  In  regions  of  low  stress. 
Because  of  cost,  however,  such  elaborate  design  is  usually  possible  only 
for  sophisticated  structures,  e.g.,  military  aircraft.  While  fastener 
problems  are  significant  In  these  Structures,  It  Is  expected  that  their 
Importance  will  increase  as  composites  usage  spreads  to  less-sophisticated 
applications. 

The  fastening  of  composites  poses  some  special  problems.  In  the  fastening 
of  metals,  normal  metal  plasticity  acts  as  an  "averaging”  mechanism 
tending  to  relieve  both  high  local  stresses  and  uneven  load  distribution 
in  fastener  patterns.  In  general,  composites  are  considerably  more 
brittle,  more  sensitive  to  stress  concentrations,  and  have  much  lower 
strain  to  failure.  Fastener  techniques  which  distribute  bearing  loads 
more  evenly  (reduce  stress  concentrations)  result  In  higher-strength 
joints.  Techniques  which  distribute  bearing  loads  Include  fastener 
arrays.  Interference-fit  fasteners,  glued  fasteners,  and  hybrid  adhesive- 
mechanical  joints.  A  typical  example  Involves  replacing  lead-aligned 
graphite  fibers  near  the  fastener  hole  with  lower-modulus  glass  fibers. 
Stress  concentration  relief  In  the  Interference-fit  fastener  case  Is 
derived  f»*om  a  matrix  "softening"  In  the  vicinity  of  the  hole.  This 
softening  Is  actually  an  Inter-or  Intra-ply  separation  which  can  occur 
without  fiber  breakage. 

1.2.  Project  Summary 

The  contractors  conducted  a  one-year  study  of  the  mechanics  of  fasteners 
In  composites.  It  should  be  noted  that  the  analytical  and  experimental 
tools  used  in  this  study  are  fairly  novel,  particularly  in  the  context  of 
composite  fastening,  and  In  that  sense  are  unique  contributions  In  their 
own  right.  The  theoretical  effort  concentrates  on  the  extension  of 
boundary  element  methods  to  determine  stresS-straln  fields  for  complex 
multiply-connected  and  three-dimensional  geometries  with  anisotropic 
materials.  Experimental  approaches  Include  a  Moire  method  involving 
high-resolution  grating  replication,  Fourier  optical  processing  of  grating 
replicas,  and  digital  data  reduction.  , 

This  report  describes  the  methodology  developed  for  both  the  analytical 
and  experimental  studies.  The  techniques  developed  were  used  to  determine 


11 


stresses  and  strains  In  the  vicinity  of  a  single-pin  loose-fit  lap  joint 
In  an  actual  glass-epoxy  composite;  these  results  are  reported.  The 
properties  of  the  composite  were  measured  so  that  theory  and  experiment 
could  be  compared.  The  comparison,  which  is  described  in  Section  4.0. 
showed  good  agreement  for  the  fundamental  case  studied. 

While  this  one-year  effort  was  successful  In  terms  of  methodology  and 
results,  it  was  viewed  from  the  start  by  the  researchers  Involved  as  the 
early  steps  In  a  more  comprehensive  research  program.  A  unique  aspect  of 
this  program  Is  that  the  experimental  and  analytical  techniques  are 
simultaneously  being  developed  and  used  by  researchers  In  the  same  group. 
This  work  Is  necessary  before  systematic  fabrication  of  structures  from 
conposltes  can  be  successful  In  engineering  and  economic  terms.  The 
effort  Is  considerably  enhanced  by  having  the  theoretical  and  experimental 
researchers  together.  Significant  progress  has  been  made  after  the  close 
of  the  one-year  contract  because  the  contractors  were  comml tted  to  the 
ideas  and  the  people  Involved.  The  researchers  feel  strongly  that  the 
major  benefits  will  come  from  research  yet  to  be  done,  and  they  suggest 
that  continued  support  Is  appropriate  for  proper  realization  of  the 
resources  and  personnel  already  developed. 

2.0.  ANALYTICAL  DEVELOPMENT,  RESULTS  AND  COMPARISON 

2.1.  Problem  Statement 

This  section  presents  the  analytical  develops*  nt  of  a  boundary  element 
formulation  for  calculating  stresses  and  deformations  In  mechanically 
fastened  composites.  Toe  purpose  of  this  development  results  from  the 
following  observations.  In  a  vast  majority  of  the  composites  literature 
In  which  stress  and/or  deformation  calculations  are  required,  the  finite 
element  method  (FEM)  Is  used.  This  method  Is,  of  course,  very  powerful 
and  has  the  advantage  of  access  to  well  developed  codes.  It  Is  not 
necessarily  the  most  efficient  numerical  procedure,  however.  Numerous 
authors  have  shown  (l)*  that  the  boundary  element  method  Is,  for  certain 
classes  of  problems,  considerably  more  efficient  than  the  FEM.  The  main 
purpose  of  this  section  is  to  develop  a  boundary  element  method  (BEM)  code 
pertinent  to  a  mechanically  fastened  composite  structure.  This  code  will 
then  be  directly  compared  with  a  current  state-of-the-art  FEM  code  to  see 
If  significant  Improvement  In  efficiency  Is  possible.  For  comparative 
purposes  both  codes  are  run  on  the  same  computer,  a  Prime  250.  In  terms 
of  problem  selection  for  comparative  purposes  there  are  a  number  of 
possible  boundary  value  problems  to  look  at  In  the  composite  fastening 
area.  ,  Examples  Include  multiple  fasteners.  Interference  fit  fasteners, 
softening  strips,  etc.  For  this  example,  a  simple  lap  Joint,  single  pin, 
loose  fit  connector  Is  considered,  and  results  using  both  numerical 
methods  are  compared.  If  Improvements  can  be  shown  In  the  simple  case.  It 
Is  reasonable  to  expect  Improvement  In  more  complicated  situations. 


♦Numbers  in  parenthesis  refer  to  references  given  at  the  end  of  Section 
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However,  code  development  for  the  more  corrpllcated  situations  are  covered 
under  future  proposed  research.  The  next  subsection  will  explain  the  BEM 
formulation.  This  is  followed  by  the  solution  of  some  example  problems 
and  a  comparison  of  BEM  and  FEM. 

2.2.  BEM  Formulation 

Let  an  anisotropic  body  (Figure  2-1)  occupy  a  finite  open  plane  region  D 
bounded  by  a  single  smooth  contour  30  which  admits  a  representation  In  the 
form  x^  *  xf(s).  The  parameter  s  is  the  length  along  3D  >rom  an  arbitrary 
origin,  and  xj  are  cartesian  coordinates  (see  Figure  2-1).  For  the 
well-known  mixed  boundary  value  problem  of  anisotropic  elastostatics: 

SaBy5  uY#$g  *  0  in  D 

on  30t  (1) 

on  30u 

30  ■  30^  +  3 Du 

Uy  represents  the  displacement  coxponents,  the  comma  denotes 
differentiation  with  respect  to  the  arguments  after  the  comma,  a$y  1*  * 
component  of  stress,  tg  is  the  corresponding  traction  vector,  nT  is  a 
component  of  the  unit  outward  normal  to  the  boundary  3D,  and  SagY$  denotes 
the  "stiffness  tensor"  for  the  material.  To  solve  this  problem  by  an 
indirect  boundary  element  method,  one  can  apply  the  Initially  unknown 
layer  of  body  forces  RY  to  the  boundary  of  the  embedded  region  and  use  the 
principle  of  superposition  to  obtain  (2): 

«ag(E)  *  /  Ha8y  (5**' )Ry(x' )ds  (2) 

Ufl(E)  ■  /  U8y  (E,x')RY(x')ds  (3) 

where  E*(Ex,E 2)  *  *'*(x*  ,y')  are  field  point  and  source  point  respectively. 
U8y  and  Hagy  are  the  fundamental  displacement  and  stress  solutions*  for  a 
point  load  in  the  y  direction  In  an  Infinite  medium.  As  £  approaches  a 
boundary  point  from  Inside  the  region,  equation  (2)  reduces  to  an  integral 
equation  of  the  second  kind  for  which  the  integral  is  defined  In  a  Cauchy 
principal  value  (CPV)  sense.  For  equation  (3)  the  corresponding  boundary 
Integral  Is  of  the  first  kind  and  need  not  be  considered  as  a  CPV.  For 
simplicity,  a  general  boundary  Integral  equation  covering  several  cases 
can  be  written  as: 

Pg(s)  ■  s  Ms>  ♦  /  g6y  (*•**)  aY(s')ds'  (♦) 

30 

8,y  -1,2 


a0YnY  *  *8 
Ug  ■  Ug 

a,8,Y  «  1,2 


are  given  in  detail  in  Appendix  A. 
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we  consider  three  cases: 

(I)  The  pure  displacement  problem 

a  3  0 

PS  3  ug 
G$y  3  Ugy 

(II)  The  pure  traction  problem 

a  »  1/2 

PS  * 

GgY  3  H0gy  n„ 

(III)  The  mixed  problem  (l.e.,  displacements  are  prescribed  on  3DU  and 
tractions  are  prescribed  on  3Dt)  requires  an  appropriate 
combination  of  case  (1)  and  case  (11). 

In  the  numerical  solution  of  equation  (4)  one  can  replace  30  by  N  straight 
line  segments  3Dj,  j  *  l,  . ...  M  on  which  the  point  loads  Ry(st)  are 
approximated  by  piecewise  constant  functions  RT(s«),  t  *  3*1.  ... 

N.  If  the  boundary  conditions  are  satisfied  at  tne  center  of  each  segment 
then  equation  (4)  reduces  to  a  system  of  2N  slmulataneous  linear  algebraic 
equations  denoted  by: 

A  R  -  b.  (5) 

where  A  Is  a  2Nx2N  coefficient  matrix,  £*(&!, R2»...R2t|)^  are  the  unknown 
point  loads  and  b3(bj,b2...b2N)T  are  th?  prescribed  boundary  conditions. 
The  essence  of  tlfe  computation  Is  the  construction  of  A.  From  the 
discretized  version  of  equation  (4)  It  can  be  seen  that  A  Is  N  blocks  of 
2*2  matrices  having  elements: 

1/2«8y  1-3 

(a  6Y)l  j  *  /  ^o$y  *  *  3  (6)  _ 

30j 

<**6y>33  1  *  3 

(agY)ij  3  /  U8  (s^s')ds'  1*3  (7) 

30  J 

where  (a*$Y),  S,Y  ■  1,2  have  been  computed  analytically  using  Figure  2-2 
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for  both  Isotropic  and  orthotropic  materials.  In  Figure  2-2,  hj  is  the 
mesh  length.  As  a+0,  the  field  point  approaches  the  boundary  point  and 
the  Integrals  over  this  mesh  length  become  singular.  The  Integrals  are 
evaluated  by  first  assuming  that  the  unknown  traction  components  R^,  Rn 
are  constant  over  the  mesh  length.  R*  Rn  can  then  be  taken  outside  the 
Integrals,  and  the  resulting  integrals  car,  then  be  evaluated  analytically 
for  a*0.  The  resulting  values  of  (a*gY)jj  are  obtained  by  taking  the 
limit  as  a-HD.  The  results  are: 

(1)  For  the  Isotropic  case 

(a*n)jj  3  2Q  hj(ln(hj/2)-l)  +  hjCOS2<*j 
<a*12)jj  *  (**2l)jj  *  hjslnajcosaj 
(a*22)jj  3  2Qhj(ln(hj/2)-l)  ♦  hjs1n2aj 

with  Q  •  (3-v)/2(l+v) ,  where  v  Is  the  Poisson's  ratio  for  material. 
(11)  For  the  orthotropic  case 

(a*ll)jj  *  -[(ciAi-C3A2)cos2aj  ♦  (C4Ai-C2A2)s1n2aj]hj(ln(hj/2)-l)/2 
(**12)jj  *  “( (ciAi+C2A2"c3A2"c4Al)cosojs‘,n®j)hj^n^hj^“1J/2 
(a*2l)jj  -  (a*i2)jj 

ta*22)jj  *  -C(ciAi-C3A2)s1n2aj  +  (C4  Ai-C2A2)s1n2<*j]hj(ln(hj/2)-l)/2 

where  c's  and  A's  have  been  defined  1w  Appendix  A. 

The  Integrals  In  equations  (6)  and  (7)  have  been  computed  .v. lyrically 
using  a  four-point  Harrls-Evans  quadrature  formula  (3).  This  quadrature 
Is  useful  when  the  Integral  Is  singular.  In  this  formulation  (the 
singular  case),  1*j  has  been  excluded,  but  In  the  adjacent  segments  (near 
singularities)  this  quadrature  Is  helpful. 

Once  the  system  (5)  Is  solved  for  £,  stresses  and  displacements  at  any 
Internal  field  point  can  be  computed  from  the  discretized  form  of 
equations  (2)  and  (3).  The  boundary  element  method  is  unable  to  predict 
stresses  and  displacements  at  boundary  points.  However,  by  excluding 
singular  points  (l.e.,  evaluating  elastic  fields  analytically)  and 
employing  the  Harrls-Evans  quadrature  for  the  rest  of  the  points,  we  can 
predict  boundary  fields  (especially  the  stress  concentration  factors  In 
the  example  problems)  quite  accurately. 

The  corresponding  computer  program,  BEM  (see  Appendix  B),  Is  based  on  the 
formulation  explained  here  and  Is  written  In  FORTRAN  77. 


2.3.  Example  Problems 

A  number  of  example  problems  were  solved  using  the  BEM  program.  All 
problems  have  a  common  geometry,  namely,  the  lower  section  of  the  simple 
lap  mechanical  joint  (see  Figure  2-3.).  Making  use  of  the  symmetry  of  the 
specimen  (for  the  principal  material  and  coordinate  axis  coincident),  the 
final  geometry  analyzed  is  shown  in  Figure  2-4.  Figure  2-5.  is  a 
schematic  of  the  boundary  element  subdivision  used.  Both  isotropic  and 
orthotropic  problems  wero  solved;  however,  only  orthotropic  results  are 
presented.  All  orthotropic  problems  were  based  on  a  particular 
glass-epoxy  composite  having  the  following  compliances: 

cn  *  5.00xl0**8  i/pSi 

ci2  *  -1.05x10-8  1/p si  (8) 

c22  *  4.76x10-7  1/psi 

C33  -  1.18x10-8  1/psl 

These  compliances  were  measured  for  the  laminate  used  in  the  experimental 
strain  analysis  phase  of  this  project  (see  Section  3-8.). 

Two  items  were  investigated  in  this  limited  numerical  study.  The  first 
Involved  hole  size,  i.e.,  the  d/w  ratio,  and  the  second  was  the  effect  of 
various  boundary  conditions.  Tables  2-1.  through  2-4.  shown  with  the 
enclosed  figures  are  self  explanatory  and  give  selected  solutions  for 
these  parameters.  Full  field  stresses  were  computed  but  only  ligament 
line  (the  line  from  the  edge  of  hole  to  the  edge  of  the  specimen)  stresses 
are  presented.  These  stresses  are  particularly  useful  for  two  reasons: 
the  maximum  stress  concentration  Is  along  this  line  (at  the  hole  edge)  and 
the  appropriate  summation  of  these  stresses  (forces)  Is  a  check  on  overall 
equilibrium.  All  computations  were  done  on  a  Cyber  750  computer.  Note 
that  in  the  following  subsection  in  which  BEM  is  compared  to  FEM^  both 
programs  were  run  on  the  Prime  250  system.  This  system  had  graphics 
capability  which  permitted  graphic  display  of  computed  data. 

2.4.  Comparison  of  BEM  and  FEM 

As  discussed  earlier,  one  of  the  main  objectives  of  this  work  is  to  make  a 
comparison  of  BEM  and  FEM  numerical  procedures.  For  this  purpose  it  was 
useful  to  run  both  programs  on  the  same  computer.  The  FEM  computer  code 
was  available  on  the  Prime  250  system.  For  this  reason,  the  BEM  code  was 
adapted  to  the  Prime  250  computer  system.  For  this  contract  FEM-BEM 
comparison  has  been  done  only  fcr  the  case  of  Isotropic  material  with  the 
simple  joint  In  tension.  This  Is  the  isotropic  version  of  the  results  In 
Table  2-1.  The  comparison  was  based  on  a  common  boundary  subdivision, 
I.e.,  In  the  BEM  code  57  boundary  elements  corresponded  to  134  quadratic 
FEM  elements,  see  Figure  2-6.  For  this  case,  the  stress  concentration 
also  Is  "known"  (4).  Results  are  summarized  In  Figure  2-6.  and  below: 


TABLE  2-1. 

Stresses  along  the  ligament  line  ab 
for  an  ortho  tropic  composite. 
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TABLE  2-2. 

Stresses  along  the  ligament  line  ab 
for  an  orthotropic  composite. 

Cosine  distribution  of  tractions  on 
the  upper  half  of  the  hole,  d/w  =  .5, 
N  -  60  (boundary  subdivisions) 
d  -  0.5 
e  =  0.5 
h  =*  1.5 
w  *  1.0 
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TABLE  2-3., 

Stresses  along  the.  ligament  line  ab 
for  an  orthotropic  composite. 

Cosine  distribution  of  tractions  cn 
the  upper  half  of  the  hole,  d/v.  -  .2, 
N  *  60  (boundary  subdivisions) 
d  •  0.2 
e  -  0.5 
h  =*  1.5 
w  «  1.0 
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TABLE  2-4. 

Stresses  along  the  ligament  line  ab 

U^mo. 

for  an  orthotropic  composite. 

Fixed  pin  case,  d/w  =  .5, 
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L.  FINITE  ELEMENT  (ANSYS) 

(  d  *  .2,  c  »,  0.5,  h  =»  1.5,  v  »  1.0  ) 


N  *  134  Quadratic  Elements 
Results:  For  d/w  «  .2  K,  -  3.1102 

X 

Total  CPU  time:  177  sec. 

2.  BOUNDARY  ELEMENT  (BEM) 

(  d  -  .2,  e  -  0.5,  h  -  1.5,  w  -  1.0  ) 


N  -  57  Constant  Elements 

Results:  For  d/w  ■  .2  Kj  ■  3.  0951 

Total  CPU  time:  27  sec. 

FIGURE  2“ 6.  FEM  and  BEM  Subdivision  Schemes  for  the  Mechanical  Joint 
Problem. 


FEM 


EXACT 


BEM 


Kj  -  3.1102  Ki  *  3.1352  Kj  »  3.0951 

CPU  time  *  177  sec.  CPU  time  *  27  sec. 

Note  that  BEM  represents  a  significant  reduction  In  computer  time  over 
FEM.  This  savings  depends  to  some  extent  on  the  particular  case,  how  many 
field  points  are  required,  and  so  on. 

3.0.  EXPERIMENTAL  METHODS  AND  RESULTS 

3.1.  Introduction 

The  experimental  components  of  this  Investigation  encompassed  two  phases. 
Techniques  appropriate  to  the  problem  and  the  material  were  developed. 

Some  studies  of  composite  pin-loaded  specimens  were  carried  out.  In 
actuality,  both  phases  ran  simultaneously,  so, the  program  was  rather 
complex  and  open-ended. 

Tha  primary  experimental  method  uses  Moire  interference  with  Fourier 
optical  processing  of  grid  photographs.*  Rugged  grids  are  applied  to  the 
specimen  by  vacuum-deposition,  and  these  grids  are  recorded  using 
high-resolution  techniques  for  each  state  of  the  specimen.  Sequential 
recording  of  grating  replicas  and  subsequent  coherent  optical  processing 
to  obtain  enhanced  fringe  patterns,  often  with  multiplied  sensitivity, 
allow  quantitative  comparisons  between  any  two  state*  of  the  specimen  at 
any  later  time.  Such  a  procedure  Is  especially  advantageous  In 
nonreversible  structural  testing. 

A  publication  (5)  discussing  the  theoretical  basis  and  showing  an 
application  of  this  Moire  technique  Is  duplicated  in  Appendix  C.  Complete 
descriptions  appear  In  Air  Force  reports  and  papers  by  Cloud’s  group  (6  to 
12)  so  further  description  Is  not  needed  here. 

It  Is  worth  noting  that  this  version  of  the  Moire  method  can  be  used  In  a 
"dirty*  environment.  It  Is  tolerant  of  poor  quality  or  course  grids  and 
poor  photographic  reproduction  of  the  gratings.  The  approach  has  been 
used  In  difficult  situations  Involving  high  temperatures  (up  to  1,5C0"F) 
long  times  (1,000  hours)  and  in  the  presence  of  convection  currents, 
vibrations,  and  so  on  (6).  The  procedures  appear  to  have  great  promise 
for  application  In  field  Investigations  of  displacement  and  strain  In 
structures  of  any  material. 

It  soon  became  clear  that  the  Moire  technique  as  summarized  above  has  a 
strain  sensitivity  which  Is  marginal  for  wort  on  composites,  as  was 
expected.  Consequently,  while  this  method  was  being  used  on  the  first 


♦In  this  report,  grating  means  a  set  of  parallel  lines,  while  grid  means  a 
pair  of  orthogonal  gratings  -25  sets  of  lines  at  90*. 
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specimens  in  this  project,  a  program  to  extend  the  sensitivity  by  a  factor 
of  10  or  more  was  begun.  Both  endeavors  were  successful.  In  order  to 
compare  theoretical  and  experimental  results  in  this  project,  the 
constitutive  properties  of  the  material  studied  had  to  be  measured.  This 
aspect  of  the  experimental  work  was  also  successful. 

This  section  describes  first  the  specimen  materials  and  preparation. 
Including  the  creation  of  Moire  grids.  The  essential  details  of  the  Moire 
process  leading  to  plots  of  strain  are  then  outlined.  Results  obtained  In 
this  way  are  given,  and  interesting  aspects  of  these  data  are  discussed. 
The  experiments  to  determine  material  properties  for  the  composite  are 
described  and  the  results  given.  Finally,  the  high-sensitivity  techniques 
are  described  In  brief,  and  some  results  which  Illustrate  the  promise  and 
flexibility  of  the  method  are  reported. 

3.2  Specimen  Fabrication  and  Material  Specification 

The  material  used  for  this  investigation  was  fiberglass-epoxy  laminate 
with  woven  fibers  (R1500/1581,  13  plies,  0.14  In.  thick)  supplied  by 
CIBA-GEIGY,  Composite  Materials  Department,  10910  Talbert  Avenue,  Fountain 
Valley,  California  92708. 

The  dimensions  of  the  specimens  are  shown  In  Figure  3-1.  It  has  been 
shown  by  Horgan  (13)  that,  when  working  with  composites,  the  end  effects 
persist  over  distances  of  the  order  of  several  widths  of  the  specimen. 

This  stress  channelling  does  not  seem  to  be  so  severe  In  this  case  based 
on  the  results  of  the  Moire  patterns  obtained,  probably  because  of  the 
reinforcing  In  the  transverse  direction.  To  be  safe,  however,  the 
specimens  were  designed  to  be  quite  long  and  narrow. 

The  specimen  was  cut  with  the  fiber  oriented  along  the  direction  of  the 
loading.  Grating  application  techniques  are  described  In  section  3.3. 
Subsequent  to  applying  the  gratings,  the  fiducial  marks,  Identlfers,  and 
code  marks  were  applied  with  Presstype  lettering.  The  specimen  was  then 
ready  for  recording  a  baseline  Moire  grating  photograph,  loading  through 
the  hole  by  means  of  a  pin,  and  recording  at-strain  grating  Images. 

3.3.  Specimen  Gratings  J _ 

Application  of  the  Moire  effect  to  any  problem  depends  on  the  successful 
deposition  of  line  grids  (or  dots)  on  the  specimen  material. 

The  photoresist  approach  to  creating  grids  on  the  specimens  was  chosen 
because  It  Is  fairly  simple,  requires  minimal  special  equipment,  and  Is 
well  proven  In  the  contractors*  laboratory.  Photoresists  for  Moire 
applications  have  oeen  studied  and  described  In  detail  by  luxmoore, 
Hollster  and  Hermann  (14,15,16).  Cloud  and  co-workers  (6  to  12)  have  used 
this  approach  successfully. 

The  photoresist  chosen  was  Shipley  AZ1350J  provided  by  the  Shipley  Co., 
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Newton,  Mass.  The  companion  thinner  and  developer  were  purchased  with  the 
resist. 

It  is  desired  for  Moire  work,  as  with  most  other  photoresist  usage,  that 
the  resist  coating  be  thin  and  uniform.  Common  application  methods 
include  spinning,  dipping,  spraying,  wiping  and  roller  coating.  The 
spinning  and  wiping  techniques  were  found  deficient  In  that  they  always 
left  some  build-up  near  the  hole  of  the  fastener  boundary,  that  Is,  In  the 
region  of  greatest  interest. 

Attention  settled,  therefore,  upon  the  spraying  method.  An  artist's 
airbrush  was  obtained  and  a  spraying  technique  which  gave  satisfactory 
uniformity  and  coating  thickness  was  worked  out  by  trial  and  error. 
Superior  results  were  obtained  by  thinning  the  resist.  Testing  was 
conducted  to  establish  a  balance  of  reslst-thlnner  proportions,  air 
pressure,  airbrush  nozzle  opening,  spraying  distance,  and  brush  motion. 

In  order  to  produce  coatings  of  the  desired  thickness,  the  photoresist 
required  thinning.  The  proportions  arrived  at  through  trial  and  error 
were  by  volume,  one  part  AZ13S0J  to  two  parts  AZ  thinner.  The  air 
pressure  provided  by  "canned  air”  sold  in  art  supply  stores  was 
satisfactory  for  the  Paasche  type  H-3  airbrush;  The  best  nozzle  setting 
for  the  airbrush  used  was  1-1/2  to  ♦  full  turns  open  from  the  closed 
position.  The  spraying  procedure  which  was  developed  called  for  laying 
the  clean  and  dry  specimen  Inclined  at  about  80"  to  the  horizontal.  The 
airbrush  containing  the  resist  was  held  about  12  in.  from  the  specimen. 
Flow  of  the  atomized  resist  was  begun  and  allowed  to  stabilize  for  about 
one  or  two  seconds,  after  which  the  spray  was  quickly  shifted  onto  the 
specimen.  At  the  range  the  air  flowed.  It  was  necessary  to  sweep  the 
brush  from  left  to  right  to  cover  the  whole  surface.  Care  was  taken  to 
assure  that  the  spray  fan  did  not  overlap  the  previously  applied  wet 
coating,  otherwise  small  bubbles  of  photoresist  started  to  form,  yielding 
a  nonuniform  surface.  It  was  absolutely  necessary  to  start  the  spray  well 
before  bringing  It  to  bear  on  the  specimen,  as  some  coarse  droplets  are 
expelled  at  the  beginning  of  flow. 

Coating  thickness  was  controlled  by  spraying  time  and  number  of  strokes 
(or  coatings).  Best  results  were  obtained  with  one  single  layer  of 
resist. 

If  the  coating  did  not  appear  satisfactory.  It  was  removed  using  either 
thinner  or  acetone.  This  operation  was  performed  as  quickly  as  possible 
to  avoid  dissolving  the  surface  of  the  specimen.  The  surface  was  then 
Immediately  washed  with  running  water  to  get  rid  of  any  excess  of  acetone. 
After  the  specimen  dried  at  room  temperature,  another  coating  of 
photoresist  was  applied.  The  coated  specimens  were  placed  Inside  a 
light-tight  container  to  await  exposure  and  development  of  the  grating 
Image. 
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3.4.  Printing  Grating  onto  Specimen 

The  Moire  grid  was  printed  on  the  photoresist  coating  on  the  specimen  by  a 
simple  contact  printing  procedure  In  which  a  grid  master  was  held  in  close 
contact  with  the  sr  :c1men  and  the  assembly  exposed  to  ultraviolet  light 
from  a  Mercury  Ian*  .  Figure  3-2.  summarizes  the  process  of  applying  the 
grids. 

The  master  grid  used  to  print  the  grid  on  the  specimen  was  made  by  using  a 
fine  metal  mesh  with  an  orthogonal  array  of  holes.  In. this  study.  Nickel 
mesh  with  2,000  lines  per  in.  was  used.  First,  a  Nickel  mesh  was  held  In 
close  contact  with  a  piece  of  flat  glass  by  spreading  soap  solution  over 
It  and  then  removing  the  surplus  with  filter  paper.  Then  the  corners  of 
the  piece  of  mesh  were  fastened  to  the  glass  with  adhesive  tape. 

The  Mercury  lamp  which  was  used  has  a  power  of  200  Matts.  The  distance 
from  the  lamp  to  the  master  grating-specimen  assembly  was  20  In.  The  time 
of  exposure  was  30  seconds.  The-newly  exposed  photoresist  was  developed 
according  to  manufacturer's  Instructions  in  the  standard  Shipley  AZ 
developer  diluted  with  water. 

Next  the  specimen  with  Its  photoresist  grid  was  placed  In  a  vacuum 
deposition  unit  (Denton  D.V.  502  high  vacuum  evaporator)  and  a  film  of 
aluminum  was  deposited  over  the  whole  surface.  The  result  Is  a  relief 
grid  in  the  metal  coating  with  excellent  reflectivity.  Aluminum  was 
chosen  because  of  Its  ability  to  resist  tarnishing.  Its  high  reflectivity 
in  thin  films,  and  its  low  cost. 

3.5.  Grid  Photography 

The  complete  state  of  strain  throughout  an  extended  field  can  be 
determined  from  Moire  fringe  photographs  obtained  through  superposition  of 
a  submaster  grating  with  deformed  and  undeformed  (baseline)  specimen  grid 
replicas.  Such  superposition  yields  baseline  Moira  fringes  and  data  (at 
strain)  fringe  patterns. 

The  set  up  used  to  accomplish  the  high  resolution  photography  of  the 
specimen  grid  Is  sketched  in  Figure  3-3.  The  camera  used  was  a  Horseman 
4X5  bellows  model.  The  lens  was  a  Carl  Zeiss  S-planar  with  focal  length 
of  120mm  and  a  maximum  aperture  of  5.6,  The  system  rested  upon  a  granite 
optical  table  and  the  camera  was  set  up  to  give  a  magnification  factor  of 
one.  The  specimen  was  placed  In  the  loading  frame  (loaded  under  tension). 
The  Tight  sources  were  two  flash  lamps  which  were  activated  by  an 
electronic, triggering  mechanism. 

Focus  of  the  specimen  Image  was  very  critical  In  this  high  resolution 
situation.  The  ground  glass  of  the  camera  was  not  satisfactory  for  this 
critical  work  because  It  was  too  coarse  and  because  such  focus  plates  are 
often  not  exactly  In  the  photoemulsion  plane.  For  focusing,  a  blank  plate 
of  the  thickness  and  type  used  In  the  photography  was  developed  and  fixed 
and  then  mounted  In  a  4X5  plate  holder  which  had  the  separator  removed. 
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SPR  s  Shipley  Photo  Resist 


Figure  3-2.  Grating  production  process. 


The  of  the  specimen  In  the  emulsion  was  examined  with  a  160X 

micros"'^  which  had  been  focused  first  on  the  image  of  the  specimen 
surface  this  process  facilitated  bringing  emulsion  and  image  Into  the 
same  ]>  t - 

After  checking  that  the  Image  of  the  specimen  was  perfectly  focused  over 
Its  entire  area  and  that  the  desired  magnification  was  correct,  the  camera 
was  locked  in  place  to  avoid  losing  the  focus  during  the  loading  of  the 
photo-plate. 

The  right  exposure  was  determined  by  trial  and  error  using  Kodak 
high-speed  holographic  film  (type  SO-253,  4X5  In.).  After  getting  the 
best  grid-replica  from  Kodak  film,  the  data  and  baseline  grids  were 
recorded  by  using  Kodak  high  speed  holographic  plates  (type  131-02,  4X5 
In.).  In  order  to  get  the  right  exposure,  the  maximum  aperture  was  used 
to  avoid  losing  sharpness  of  the  grating,  and  the  Intensity  of  the  flash 
lamps  was  reduced  by  adding  ground  glass  pieces  in  front  of  the 
flashlamps.  These  glasses  cut  down  the  Intensity  of  the  light  by  roughly 
one  f-stop  per  piece.  Also*  It  Is  worth  noting  here  that  the  angle  of 
incidence  of  the  Illumination  was  chosen  by  trial  and  error  to  give  the 
best  contrast  In  the  grating  Image.  Shadows  of  the  three-dimensional  grid 
structure  evidently  play  an  Important  role  In  grating  visibility.  The 
exposed  131-02  plates  were  Individually  developed  In  Kodak  HRP  developer 
for  three  minutes  and  fixed  In  Kodak  fixer  for  the  same  amount  of  time. 

After  completion,  each  grating  plate  was  examined  and  Inspected  for 
diffraction  efficiency  to  assure  that  the  photography  had  been  successful. 
It  was  labelled  with  specimen  number,  loading,  and  magnification 
conditions  for  future  reference  and  stored  In  a  rack  to  await  optical 
construction  of  Moire  fringe  patterns. 

3.6.  Summary  of  Optical  Processing  to  Create  Fringe  Patterns 

The  grating  photographs  of  this  experiment  produced  an  assembly  of 
photographic  plates  of  the  undeformed  (baseline)  and  deformed  (data) 
specimen  gratings  as  well  as  a  submaster  grating  having  an  Integer 
multiple  of  the  specimen's  grating  spatial  frequency.  The  creation  of 
Moire  fringe  patterns  from  these  plates  and  the  reduction  of  Moire  fringe 
data  have  been  described  In  detail  (9  to  12).  The  steps  required  to 
produce  Moire  fringe  photographs  from  these  grid  records  were  as  follows: 

1.  A  photoplate  of  the  undeformed  specimen  grid  was  superimposed 
with  a  master  grating  having  a  spatial  frequency  of  1000  Ipl 
(which  was  half  the  spatial  frequency  of  the  specimen  grating) 
plus  or  minus  a  small  frequency  mismatch. 

2.  The  superimposed  gratings  were  clamped  together  and  placed  In  a 
coherent  optical  processor  and  adjusted  to  produce  a  correct  base 
line  (zero  strain)  fringe  pattern  at  the  processor  output,  where 
ft  was  photographed.  This  photograph  corresponded  to  the 
superposition  of  the  2nd  diffraction  order  of  the  master  grating 
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with  the  first  diffraction  order  of  the  corresponding  specimen 
grating. 

3.  Steps  1  and  2  were  repeated  for  the  other  component  of  the 
,  specimen  grid  l.e.  the  grating  at  90*  to  the  first  one. 

4.  Steps  1  to  3  were  repeated  with  the  photographs  of  the  deformed 
grid  in  order  to  create  the  "data"  or  "at  strain"  fringe 
patterns.  The  same  submaster  plate  was  used. 

5.  The  fringe  patterns  were  enlarged  and  printed  with  high  contrast 
In  a  convenient  size  equivalent  to  about  5  times  the  actual 
specimen  dimensions. 


6.  The  prints  were  sorted  and  coded  for  Identification. 

7.  Computer  digitizing,  data  reduction,  and  plotting  were  performed 
on  each  photograph. 

8.  Digital  processing  of  Moire  data  finished  the  analysts. 

Digitization  of  the  fringe  patterns  and  subsequent  data  reduction  followed 
the  procedures  described  in  detail  by  Cloud  and  colleagues  (9  to  12). 

Some  of  the  essentials  are  reproduced  in  Appendix  C. 


3.7  Experimental  Results 


In  this  Investigation,  the  specimen  had  a  two-way  grating  (a  grid).  By 
using  the  optical  processor,  a.  Moire  fringe  pattern  was  formed  separately 
for  each  direction.  One  pattern  was  formed  to  get  horizontal  fringes 
(perpendicular  to  the  direction  of  loading),  and  this  fringe  pattern  was 
used  to  measure  strain  In  the  direction  of  loading.  The  other  pattern  was 
formed  with  vertical  fringes  (same  as  direction  of  loading);  fringes  In 
this  direction  were  used  to  measure  transverse  strain  (perpendicular  to 
the  direction  of  the  load).  The  photograph  of  the  fringe  pattern  was 
recorded  separately  for  each  direction  and  for  each  loading  step.  Sample 
photographs  of  the  Moire  fringe  patterns  obtained  from  the  composite 
material  fastener  specimen  for  the  horizontal  and  vertical  directions  are 
shown  In  Figures  3-4.  to  3-7. 


Figures  3-8.  through  3-11.  show  the  measured  strain  ex  and  ey  for  several 
different  lines  on  the  specimen.  The  whole-field  nature  of  the  Moire 
method  means  that  a  great  deal  of  data  are  generated.  Such  data  are 
always  difficult  to  present  and  assimilate.  The  plots  shown  here  are  a 
reasonable  compromise  between  completeness  and  confusion.  The  next  step 
Is  to  develop  strain  contour  maps,  wnlch  would  be  easier  to  understand  at 
a  glance.  Such  a  step  was  not  within  the  scope  of  this  project.  Programs 
for  three  dimensional  computer  graphics  representations  of  the  strain 
contours  are  being  developed  by  the  Investigators. 
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Figure  3-6.  Photograph  of  Moire  pattern  showing  displacements 
parallel  to  load  lint  for  /ero.  load  with  pitch  mismatch. 
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ligaments  of  fastener  hole. 


Figure  3-9.  Strain  perpendicular  load  to  axis  (ex)  along  several  lines  In 
bearing  region  near  hole. 
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3.8.  Experimental  Evaluation  of  Elastic  Constants  for  the  Composite 

In  order  to  compare  theory  and  experiment*  It  Is  necessary  to  know  with 
precision  the  material  properties  of  the  composite  usede  An  experiment  to 
determine  these  parameters  was  conducted. 

There  are  two  goals  In  the  design  and  test  of  a  tensile  specimen.  First, 
the  existence  of  a  statically  determinant,  uniaxial  state  of  stress  within 
the  test  section  must  be  assured;  however ^  producing  such  a  state  of 
stress  In  the  laboratory  Is  not  a  trivial  task.  Several  analytical 
studies  have  revealed  that  the  first  goal  may  be  accomplished  by 
establishing  the  specimen  geometry  such  that  the  length-to-wldth  ratio  Is 
a  practical  maximum  (13,  17,  18).  In  the  current  study  the  maximum  length 
was  chosen  to  be  12  In.  and  the  width  one  In.  The  second  goal  of  the 
design  and  test  of  the  tensile  specimen  Is  to  assure  that  the  elastic 
responses  In  the  In-plan?  shear  and  traverse  tension  modes  are  constant 
and  that  failure  will  occur  within  the  specimen  test  section. 

Determination  of  Ej,  £2  and  v  Is  quite  straight  forward.  Two  tensile 
tests  are  required,  and  for  this  purpose  two  specimens  were  constructed: 
one  with  the  warp  fibers  oriented  along  the  direction  of  the  loading  (for 
determination  of  Ex  and  vx2)  and  the  other  with  woven  (weft)  fibers  In  the 
direction  of  loading  (for  determination  of  Eg  and  V21). 

As  shown  In  Figure  3-12.,  four  single-element  strain  gages  of  type 
EA-13-075AA-120  from  Ml cromeasurements  were  used  to  check  the  uniformity 
of  the  stress  field.  Strain  gage  rosettes  on  both  faces  ef  the  specimen  . 
were  also  used  (type  CEA-06-062UR-120) ,  first  to  determine  any  effect  of 
bending  on  the  specimen  and  then  to  obtain  the  measurements  needed  to 
perform  the  calculations.  Effects  of  transverse  sensitivity  were  checked 
and  found  negligible. 

The  specimens  were  mounted  In  a  loading  frame  using  grips  especially 
designed  to  avoid  any  clamping  of  the  ends  and  to  allow  rotations,  thus 
avoiding  any  end  effects. 

Every  specimen  was  loaded  In  Increments  up  to  500  lbs  and  then  unloaded 
and  reloaded  to  the  same  stress.  The  tests  were  at  room  temperature 
(75*F) .  The  strain  readings  were  made  with  a  digital  strain  Indicator 
from  Northern  Technical  Services,  Inc. 

Figure  3-13.  shows  typical  plots  of  stress  vs  strain  obtained  from  this 
experiment. 

The  results  obtained  for  Young's  modulus  and  Poisson's  ratio  are  as 
follows: 

Ex  -  3.188  x  10*  psl  E2  ■  3.0824  x  10*  psl 

vj2  •  0.11  V2x  ■  0.11 
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Figure  3-12.  Specimen  showing  location  of  strain  gages. 
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A  test  which  Is  planned  for  subsequent  study  is  to  use  an  Instron  machine 
to  perform  the  tensile  test  under  different  room  humidity  and  room 
temperatures.  It  Is  suspected  that  humidity  and  temperature  affect  the 
material  properties  of  this  specific  composite. 

3.9  Experimental  Results  and  Discussion 

A  typical  composite  pin-loaded  joint  exhibits  a  strong  coupling  between 
axial  strain  and  bearing  capacity  at  any  given  pin  location. 

In  this  Investigation,  the  specimen  had  a  two-way  grating.  By  using  the 
optical  processor,  the  Moire  fringe  patterns  were  formed  separately  for 
each  direction.  One  pair  (Figures  3-4.  and  3-5.)  was  formed  with  vertical 
fringes  (direction  parallel  to  the  load)  and  these  fringe  patterns  were 
used  to  measure  strain  ex  In  the  traverse  direction  (perpendicular  to  the 
direction  of  the  load).  The  other  pair  (Figures  3-6.  and  3-7.)  were 
formed  with  horizontal  fringes  (perpendicular  to  the  loading  direction); 
fringes  In  this  direction  were  used  to  measure  strain  ey  in  the  same 
direction  as  the  load. 

Both  the  photographs  of  the  horizontally-oriented  (from  the  horizontal 
grating)  fringe  patterns  and  the  vertically-oriented  (from  the  vertical 
grating)  fringe  patterns  yielded  some  Interesting  results. 

I.i  Figure  3-8.  and  Figure  3-9.,  plots  of  ex  vs  position  are  shown.  In 
Figure  3-10.  and  Figure  3-11.  ey  vs  position  are  shown. 

From  Figures  3-10.  and  3-11.  we  can  notice  that.  In  the  area  of  contact 
between  pin  and  composite,  ey  Is  highly  compressive;  but,  as  we  move  away 
from  the  edge  of  the  hole  towards  the  end  of  the  specimen.  Its  value 
decreases  to  an  average  value. 

Now,  looking  at  the  right  hand  edge  of  the  hole,  is  of  tensile  nature; 
as  we  move  away  from  the  edge  of  the  hole  towards  the  edge  of  the  specimen 
It  decreases  to  an  average  value  (shown  by  lines  x3  to  x7  In  Figure 
311.). 

Figure  3-14.  shows  a  distribution  of  strain  based  on  the  results  from  the 
the  Moire  fringes.  This  figure  Is  developed  from  Figures  3-10.  and  3-11. 
plus  similar  strain  plots  created  for  additional  lines  on  the  specimen. 

Much  more  Information  useful  to  the  designer  can  be  extracted  from  these 
Moire  photographs.  As  an  example,  examine  the  areas  shown  In  Figure  3-14. 
where  the  normal  strain,  l.e.,  ey,  changes  from  tensile  to  compressive. 
Along  that  interface,  we  find  a  very  high  shear  strain  which  In  some  cases 
cur  produce  delamination  and/or  failure  by  shearing  *  ‘especially  along  the 
transition  zone  from  high  compressive  strain  to  tensile  strain  near  the 
fastener.  Such  features  of  the  Moire  results  allow  the  designer  to  get 
very  useful  Information  anywhere  on  the  surface  of  the  specimen  providing 
better  criteria  for  the  determination  of  the  optimum  combination  of 


parameters  such  as  distance  from  hole  to  edge,  ratio  of  pin  to  hole 
diameters,  fiber  orientation,  lay-up  sequence,  etc. 

3.10.  High-sensitivity  Interferometric  Moire  Technique 

From  the  results  for  ex  (the  smaller  strain  component),  it  clear  that  the 
Moire  technique  has  a  marginally  low  strain  sensitivity  for  work  on 
composite  materials.  This  Idea  Is  suggested  on  Figure  3-8.  which  shows 
poor  agreement  between  values  of  ex  at  areas  located  at  the  right  and  left 
of  the  specimen  (shown  by  lines  xO  to  x2  and  xA  to  xC  respectively). 

The  Moire  method  yields  full-field  information  of  the  In-plane  surface 
displacements.  It  has  great  potential  for  the  macroscopic  strain  analysis 
of  composites,  and  this  method  does  not  suffer  any  limitation  due  to 
anisotropy.  Inhomogeneity,  or  Inelasticity  of  composite  material. 
Successful  Moire  strain  analysis  requires  that  the  sensitivity,  which  Is 
governed  by  the  grating  frequency,  be  matched  to  a  degree  with  the 
magnitude  of  the  deformation  which  Is  to  be  measured.  Thus  an  extensive 
program  to  extend  the  sensitivity  by  a  factor  of  10  or  more  was  begun. 

An  Interferometric  technique  similar  to  that  described  by  Post  (19,  20) 
and  Walker,  Mckelvle  and  McOonach  (21)  was  developed  in  this  laboratory. 

It  Is  evident  that  this  technqlue  for  measuring  the  strain  distribution 
should  operate  over  a  sufficiently  short  gage  length  to  elucidate  the 
details  of  the  strain  distribution,  while  at  the  same  time  giving  an 
overall  picture  of  the  material  behavior.  In  this  interferometric  Moire 
technique,  the  specimen  grating  Is  a  phase-type  grating.  The  analysis  Is 
carried  out  by  using  overlapping  beams  of  coherent  light  to  create  the 
specimen  grating  and  the  master  grating  which  Is  projected  onto  the 
deformed  specimen. 

The  particular  virtues  of  the  system  developed  for  this  Investigation 
are: 

1.  The  efficiency  of  light  use  Is  very  high. 

2.  No  rigid  connection  Is  required  between  the  specimen  and  the 
system.  This  last  feature  Is  relevant  In  view  of  the  convenience 
of  performing  measurements  In  environments  which  are  not  so  Ideal 
as  a  vibration-free  optics  laboratory. 

3.  Measurements  can  be  performed  in  three  different  directions, 
yielding  a  map  of  strains  In  the  same  number  of  directions  and 
allowing  calculations  of  maximum  strains.  In  composites.  It  does 
not  make  much  sense  to  talk  about  principal  directions,  since 
they  depend  on  the  material  directions;  stilt,  information  In 
three  directions  will  provide  very  useful  Information  for 
specific  situations  of  fastener  design. 

Preliminary  testing  Of  this  technique  has  been  conducted  and  the  results 
appear  to  be  very  promising.  Figures  3-15.  and  3-16.  show  Moire  fringes 


showing  displacements  perpendicular  to  the  direction 
of  loading  for  no  load  with  rotational  mismatch. 
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of  displacement  In  the  x-di recti on  (perpendicular  to  the  direction  of  the 
loading).  Notice  the  difference  between  these  fringe  patterns  and  those 
obtained  by  the  traditional  Moire  technique  (Figures  3-4.  and  3-5.). 

4.0.  COMPARISON  BETWEEN  EXPERIMENTAL  AND  NUMERICAL  RESULTS 

Section  3  describes  the  methodology  of  obtaining  strain  in  the  actual 
composite  material  by  Moire  techniques  which  are  self-rail bra ting  and 
which  do  not  depend  on  assumptions  of  material  behavior.  Typical  results 
are  reported.  Section  2  described  a  boundary-element  approach  to  the 
problem  and  gives  results  from  that  method.  Some  finite  element  results 
are  also  computed  and  used  for  comparison  of  the  two  numerical 
approaches. 

It  Is  Instructive  and  valuable  to  designer,  theoretician,  and 
experimentalist  to  compare  the  results  from  these  widely  diverse 
approaches  to  the  problem.  Since  the  methods  are  so  different,  good 
agreement  of  results  would  imply  that  the  results  are  probably  correct. 
Disagreement  would  give  an  Idea  of  the  magnitudes  of  errors  In  either  or 
both  sets  of  results. 

To  perform  the  comparison,  some  strains  measured  by  Moire  were  converted 
to  stress  by  use  of  generalized  Hooke's  law.  This  procedure  Inserts  ah 
assumption  of  material  behavior  Into  the  experimental  findings,  but  It 
does  facilitate  comparison.  It  Is  probably  more  reasonable,  but  also  more 
difficult,  to  iise  Hooke's  law  to  convert  the  numerical  results  to  strain. 
The  conversion  of  experimental  data  was  carried  out  for  ey  at  selected 
points  along  the  ligament  line  between  hole  boundary  and  specimen  edge. 
These  stresses  were  normalized,  then  plotted.  They  are  compared  with  the 
matching  finite  element  values.  The  comparison  Is  shown  graphically  In 
Figure  4-1.  The  agreement  between  numerical  and  experimental  findings  is 
excellent,  especially  when  one  considers  the  fact  that  these  composites  do 
show  some  non homogeneity ,  l.e.,  properties  vary  slightly  from  point  to 
point. 

This  comparison  Is  not  extensive,  but  It  does  Indicate  that  the  work  and 
the  results  are  probably  correct.  Comparisons  for  other  critical  areas  of 
the  stress  field  are  planned.  In  addition.  It  would  be  advisable  to 
obtain  a  measure  of  the  Inhomogeneity  of  the  material  so  that  a  "scatter 
band"  of  expected  stress  for  a  given  load  situation  can  be  established. 

The  results  presented  suggest  that  numerical  and  experiment  results  would 
both  He  within  any  such  scatter  band. 
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L  APPENDIX  A 

• 

In 

J 

•  FUNDAMENTAL  SOLUTIONS 


V 


' 


l 

F 


V 


For  the  Isotropic  case 


The  stress  and  displacement  fields  at  a  point  (x,y)  due  to  a  unit  point 
load  acting  at  the  origin  of  coordinates  in  the  x-dlrection  in  an  Infinite 
plate  with  material  properties  G  and  v  under  plane  stress  are: 


4G  Ux 

Q  Inr2  +  y2/r2 

4G  Uy 

.  1+v 

-xy/r2 

«xx 

x(P  ♦  x2/r2)/r2 

Hxy 

y(P  +  x2/r2)/r2 

Hyy 

x(-P  +  y2/r2)/r2 

where  P  *  (l-v)/2(l+v) ,  Q  ■ 

.5  +  2P  and  r2  -  x2  +■ 

the  y-dl recti on,  the  results  are: 


4G  Ux 

-xy/r2 

4G  Uy 

Q  lnr2  +  x2/r2 

«xx 

* 

TT 

y(-P  ♦  x2/r2)/r2 

Hxy 

x(P  +  y2/r2)/r2 

Hyy 

y{P  ♦  y zlrz)/rz 

For  a  unit  load  In 


The  tractions  are: 


Tx,k  *  Hxx,k  cosa  +  Hxy,k  slna 

Ty,k  *  Hxy,k  cos<*  +  Hyy,k  s*nct  i 

where  k  «  either  x  or  y  and  a  -  angle  tha*: 
the  x-axIs. 

For  the  orthotropic  case 
For  a  unit  load  In  the  x-dl recti on 


the  exterior  normal  makes  with 


Uxx 

(cxAi  In  r2  -  C3A2  In  ri)/2 

Uyx 

AiA2(4i~42)/C22 

Hxxx 

(cir2  -  C3ri)x 

-  k 

Hxyx 

(cir2-C3ri)y 

Hyyx 

( 5lc3rl  ”  <s2clr2^ x 

For  a  unit  load  In 

the  y-dl recti on 

Uxy 

AiA2Ui-*2)/C22 

Uyy 

(C4A1  lnri  -  c2A2  lnr2>/2 

Hx*y  •  k 

(C2r2  -  C4ri)y 

Hxyy 

t5lc4rl  "  52c2r2)* 

Hyyy 

(4lC4ri  -  $2c2r2)y 

with 


k  -  l/2ir(52“«l), 

H  ■  l/(4jx24y2),  1*1,2 


Ai  -  cj2  *  «1  C22,  1-1,2 

♦l  -  *2  ■  s1n_1Cyx/rir2  (/SJ  "  *^3. 
ci  *  /SjjT  (al”Ci2/C22) » 
c2  ■  (l”^ici2/cn) , 

C3  -  /^T  {«2-ci2/c22). 

C4  m  /7i  (1-<2C12/ch). 

$1  and  «2  *r®  the  roots  of  the  characteristic  equation  of  the  material; 

C22  52  **•  (2012^33)^  +  cn  -  0 

Mote  that  i's  are  either  real  or  complex.  If  they  are  real,  then  since 

51*2  *  CII/C12  >  0,  they  are  either  both  positive  or  both  negative.  For 
clearness,  the  4‘s  are  taken  to  be  real.  For  complex  $'s,  «2 
necessarily  the  conjugate  of  61  and  the  c's  obtained  above  are  substituted 
In  (11)  to  obtain  the  equivalent  complex  forms  of  the  fundamental 
solution,  the  real  parts  of  which  represent  U‘s  and  H's. 
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APPENDIX  B  i 

THE  COMPUTER  PROGRAM  "BEM"  JJ 

s 

p" 

* 

i 


Major  steps  In  constructing  BEM  computer  program 

Step  1.  Provide  the  following  data: 

Boundary  points  coordinator  (XB,  YB) 

The  angles  that  outward  normal  at  boundary  nodes  make  with 
+x-ax1 s  (T) 

Field  points  coordinates  (XF.YF) 

Boundary  conditions  In  X  and  Y  (BC), 

BC  *  0  displacement  prescribed 
BC  »  1  stress  prescribed 

Boundary  values  (B) 

Integration  nodes  and  weights  (Z  ,W) 

Compliances  (c's) 

Step  2.  Construct  the  Influence  matrix,  A 

Step  3.  Solve  AX  *  B 

(stores  results  In  8) 

Step  4.  Calculate  displacements  and  stresses  at  given  field  points 
Remark: 

The  computer  program  BEM  Is  developed  for  real  and  &2r  11,6  complex 
version  of  BEM  follows  directly  from  the  discussion  given  In  Appendix  I. 


PROGRAM  B EM  C INPUT*  GUTPUT  > 

C  ************************  »************************************i 

PARAMET£RCN=60»N2=2*N*M=26> 

REAL  X  B (h ♦  1 )  *YB(N*1)*TCN)  *X F CM > *YF CM) *BCC N2 >  » A CN2 *N2 > *B CN2 > * 
♦ZC4),UC4> 

READ**  C1I*C12*C22*C33 

READ  *  * C  Z  ( I  ) *1=1*4)* CUCI)*I=1*4) 

R£AD**CXBCI)*YB(I>*I=1*N)*CTCI)*I=1*N) 

READ**  CBC  ( I) *  I=1*N2)« (B(I)*I  =  1«N2) 

READ** CXFCI)*YFCI)*I=1*M) 

PI=3. 141592653589 

C  COMPUTE  CONSTANTS  IN  THE  FUNDAMENTAL  SOLUTIONS 

P1=CC12*.5*C33)/C22 
-  r-  vrr--  Xl=Pl-SQRTtPl*Pl-Cll/C22> 

^  ^■-’w  :tn(3=SQR'rC)f2)-SGRTCXll  -“5:  *■.  r;  \ 

CN=2.*PI*CX2-X1) 

-w*:*-  A1=C12-X1*C22^ . — -  •?■•••  .» 

A2=C12-X2*C22  - 

C1=CX1-C12/C22J*SQRTCX2) 

C2=(1.-X1*C12/C11)*SQRT(X2> 

•r-  . C3=C X2-C12/C22 )*SQRTCX1 )  -  - -  * 

C4=C1.-X2*C12/C11)*SQRTCX1> 

C  '  CONSTRUCT  THE  INFLUENCE  MATRIX*  A 
•  -  -  00  100  I  =  1*N  - 

11=2*1-1 

■'  >L»ire.., '.•*■..•<12^2*  _T':. ;  _  „.-J  ::  .  "-»•  -  .  ......  .  ..  a, !?;.  '■■■:•'  . 

s"  CT=COSfTCir>  •  -  • 

ST=SINCTCI>> 

FX=.5*CXBCI*1)*XBCI)) 

.  *  .  Vi  FY=.a*tYBCI*l>*YB<I)) 

•  •  00  100  J=1,N  ‘  - 

'  AX=XBC J*1>-XBCJ> 

AY=Y8C J*1)-YBC J) 

H=SQRT(AX*AX*AY*AY) 

QI=-.5*F*CAL0GC*5*H>-1.) 

IF  CI.EQ.J)  GOTO  50 

UXX=UXY=UYY=HXXX=HXYX=HYYX=HXXY=HXYY=HYYY=0. 

DO  7  K=1 .4 

tv;-.'*iX=FX-AXv*z<K)-XB{j>. 

Y=FY-AY*ZCK)-YBCJ> 

R1=1./CX1*X*X*Y*Y> 

R2=1./CX2*X*X*Y*Y> 

UXX=UXX*.5*CC1*A1* ALCG(R2)-C3*A2*AL0G(R1) >*H*UCK> 

UXY=UXY +A1* A2*ASIN C  X*Y*X3*SGRTCR1*R2>  )/C22*H>*U  CK) 

•_*r;  •:  HXXX=ttXXX*<Cl*R2-C3*Rl)*X*H*ttCK) 

•;r^/ArJi>fe*>.v,vHX.YX=HX.YX45<Cl-*R2-C3*Rl)*Y*H*<y(K):" 

HYYX  =  HY YX  +  CX1 *C3*R 1-X2  *C1*R2 )*X*H*UCK) 

HXXY=HXXY*CC2*R2-C4*R1 >*Y*H*U CK) 
HXYY=HXYY*CX1*C4*K1-X2*C2*R2)*X*H*«CK) 

7  HYYY=HYYY*CX1*C4*R1-X2*C2*R21*Y*H*UCK> 

IFCBCC  ID  •EQ.O.)  GOTO  10 

AC Il*2*d-1)=HXXX*CT*CT4HYYX*ST*ST*2**HXYX*CT*ST 
AC tl,2*U)=HXXY*CT*CT*HYYY*ST*ST*2.*HXYY*CT*ST  / 

GOTO  20 

10  AC Il,2*J-l)=UXX*CT*UXY*ST 
ACI1*2*J)=LXY*CT*UYY*ST 
20  IFCBCCI2)*EQ.0.)  GOTO  30 

AC  12*2*0-1 )=CHYYX-HXXX)*CT*ST*HXYX*CCT*CT-ST*ST> 

ACI2*2*vl)  =  C  HY  Y  Y-HXXY  >  *CT*ST*NX  Y  Y*  CCT*CT-ST*ST) 

GOTO  100 

30  ACI2*2*J-1)=UXY*CT-UXX*ST 
AC  12*2* J)=UYY*CT-UXY*ST 
GOTO  100 

50  ACIl*2*d-l)=*5*CN*CBCC I 1 > > *C7*C 1 ,-BC C I 1 > > *QI *C Cl* A 1-C3* A2 > *C1 
AC  II ,2*J)  =  .5*CN*BCC  I1) *ST* C 1 .-BC C 1 1 ) ) *Q I  * CC 1  * A1 -C3 * A2 > 4ST 
<  AC  I2*2*J-1)=-.5*CN*BCC I2)*ST*C1.-BCCI2>)*QI*  CC2*A2-C4*Al)*ST 
ACI2*2*U)  =  .5*CV*CBCCI2)  >  *CT*<«  1  .-BC  C 12  >  )  * GI *  C  C4*  A1-C2*  AS  )*CT 
lOu  CONTINUE 
C 

C  SOLVE  A • X  =  B 

CALL  HaTX  <:«?.  A  *fc) 

ORINT  5oC 


v.  .  •  .  *.v 


»•  .V.s  . 

‘  ‘>V-\ 


-'-'M-v-.'Iv'c.v 
"s  '■*' 


C  CALCULATE  DISPLACEMENTS  AND  STRESSES  AT  GIVEN  FIELD  POINTS 
DO  300  1  =  1  * M 
UX=UY=TXX=TXY=TYY=0 . 

CT=CCS  C  T  < I >  ) 

ST=SIN(T(I>> 

00  *00  J=lt N 
K=2*  J 

Ax=xe<  J*»1)-XBC  J) 

AY  =  YB(  J-*1)~YB<  J) 

H=SQRT(AX*AX+AY*AY> 

IF<(XFCI>-„5*<XB<J*1)*XB<J>  >  >  *  •  2*  <  YF  <  I  > -.5*  t  YB  t  J-*-l  > -*>Y  B  C  J)  >  >**2 
♦LE.l.E-8)  GOTO  150 
00  110  L=J.*4 

X=XF<T)-AX*Z<L>-XBCv»>  - 

Y=YFCI>-ZiL)*AY-Y8<J> 

R1=1*/<X1*X*X+Y*Y> 

R2=1./(X2*X*X*Y*Y) 

UXX=.5*(C1*A1*AL0G<R2>  -C3*A2  ^LOG(Rl)  ) 
UXY=A1*A2*ASIN<X*Y*X3*$GRT(R1*R2> )/C22 
I*YY  =  .5*<C4*A1*AL0G<R1> -C2*A2 *ALOG <R2>  ) 

HXXX=(C1*R2-C3*R1)*X 

HXYX  =  CCl*R2-C3*RlHrY 

HYYX  =  ( X1*C3*R1-X2*C1*R2 )  *X 

HXXY=CC2*R2-C4*R1J*Y 

HXYY=CX1*C4*R1-X2*C2*R2>*X 

HYYY=*X1*C4*R1-X2*C2*R2) *Y 

UX=UX*(B<K-l)*UXX*B<K)*UXY)*fr*UCL) 

UY=UY*{B(K-1)*UXY*B(K) *LY  Y>  *  »■*  W  <  L> 
TXX=TXX4(B(K-1)*HXXX+B(K)*HXIY)*H*W{L> 
TXY=TXY*<B(K-1)*HXYX4B(K)*HX1Y)*H*«(L> 
TYY=TYY4(BCK-1)*HYYX+B<X)*HYYY)4H*UCLI 
lie  CONTINUE 
GOTO  200 

150  0I=-H*<ALCG<H)-1.) 

LXX=GI *<  <C1*A1-C3*A2)*CT**2-*-  { C4*  A  1-C2* A2 )  *S T**2  > 
UXY=GI*<C1*A1+C2*A2-C3*A2-C4«A1>*CT*ST 
LYY=GI*< (C1*A1-C3«A2)*ST**24 <C4*A1-C2*A2>*CT**2 > 
LX=UX4UXX*B(K-1|4UXY*B(K) 

UY=UY4UXY*BCK-1)4UYY*B(K) 

TXX=TXX*.5*CN*CCT*B<K-1)-ST*B<K> ) 
TXY=TXY4,5*CN*C-CT*B(K-1)4ST*B(K) ) 
TYY=TYY4«5*CN*CST*B<K”1)4CT*B(K>) 

,  2 C 0  CONTINUE 

3 C 0  PRINT  4  G  C  »  1  *  X F ( I ) » YF ( I ) * UX « L Y ♦ TX X ♦ TX Y « TY Y 

400  FORMAT  ( 1H 0 ,5X , 12  * 3X *2 < 4X fF8 .4 > »6X,2< 4X ,F12.8>  *6X * 3t4X , F 12 .8 >  > 
500  F0R*AT<*&*,5X,«N0*,T21«»X*,T32t*Y*»T52**UX**T68f*UY* 
♦,T9C**TXX**T106.»TXY* , fl22»*TYY*/) 
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kola  aro  than  adddi  haro  baaa  awdalaiad  by 
ai  a  aiaala  aqaoa  taqaaacy  addek  any  ba  (ha 
amino  fraqouey  or  oaa  of  ha  hanaoafaa. 
and  addin  aeabd  frtaon  wkkk  eearerooad  la 


A  dhtiecdro  famaa  of  tka  approach  ia  that  tka  output  Haa  aad  dal  prodaii  af  MOB  W  (J*  Haaafron)  tptet 

knot  of  tka  trodaoi  a  eoaddarad  id  coaaal  of  a  drababk  fwoaaary  aan  tka  tpnaro  roriaea  akh-SMUay  photo- 

u«aai  plaa  a  (real  deal  of  oihar  jafonaadoa.  Tha  da-  min.  Tha  n^roaa  proda*  «aa  pkototrepkat  ttdao 

ponaai  doaai  a  made  <wMa  by  «fU*«  it  taro  ad  the  *i«h-mohdloa  ndadqaa  before  aad  aAar  tka  odd  work 

can  iaforamioo.  Oaa  haa  a  eartaia  latitada  ia  nlacihn  "•*  •admad.  Throe  (redag  pkocoa  warn  (haa  »P*r- 


tha  nforamdoa  that  ia  non  oaafai. 

Rational!  ration  of  Torn  Approaches 
Hoeing  taro  aaplaaaaory  roodak  af  oaa  procero  rakaa 


layad  n  wad  wad  aiahar  free— ry  nihoiaatar  (rampa 
ia  aa  nytka  yreeacwr  ia  aTOar  to  obiaia  ward  frinpaa. 

Tha  qierharo  tmiaa  phaaocraehe  had  a  tpetUl  ta> 
oaaacy  of  74*  Ipi  (JO  garoriroaj  wkkk  rendu  with  a 
■paamaa  (naan  of  1009  Ipi  najaiflod  IJ  thaaa.  Thaaa 
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pattern*  acre  made  l.x  ,hertwg  purpino  amh  wtmu >mi 
Of  741  Ipt.  Fur  several  of  the  specimens,  all  or  most  of 
these  several  moirt!  patients  mere  attained  to  (am  re- 
duodaocy  of  data.  This  earn  information  eras  useful  • 
pana|  an  apprenalioa  of  pwhoMt  errors  and  ia  stud', 
the  sense  tvuy  of  the  daw  oaahty  to  rsriattoas  «  ise 

li  was  also  poanble.  at  this  suae,  to  seise,  nuaasuts 
which  had  dcasay  and  difTraaioa  charaetesieta  which 
bctanced  whh  the  properties  of  the  yoats  (rating 
raphes  to  produce  the  beat  friage  patterns.  Also,  the  ray 


gMIpttagwa  twtlrh  Iry  aaihlplltahoa  af  3.  or  of  IM» 
Ipi  foe  a  ■nhfpgraUna  of  1  Aa  adraaead  tadmique, 
whldl  is  ■satis  aid  below,  gave  a  awdUphcatton  of  4. 
The  various  patch  adataaachat  ware  chonm  ta  yield  the 
doeaet  Mats  waring  nhtelaable  whh  good  Mage  visibility. 
Maw  grating  phosoplatea  wara  procaeaad  whh  at  least 
thraa  aauaaaaeh  levels,  and  snsnatimsi  whh  atom  than  one 


tMMvky  mMUokcatiao  ftttftf  for  ckmkiot  pupom  ood 
bacanae  h  was  not  posable  In  eases  the  qnalhy  of  a  donee 
fringe  patient  through  the  camera  viewfinder,  which  was 
used  without  a  magnifier. 

To  be  spsdflc,  for  moat  of  this  study,  each  baaahae 
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222S  and  2?Jd  Ipl.  Far  specimen  gratings  of  posr  quaUty, 
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tutmmtm  of  P33S  aod  \m  ipa  Mid.  Som  fHo* 


TW  15*om  i 


i  of  <k«  frio 


MotyM.  Tfct  ratio  of  prtoird 

(0 oi  Mf-v  — )  MM  aprsoM 
pint  ike  award  no 


aa  H§.  IQL 


Opdcd  PraMn  to  taoproo  QraHoQ  PhotOQnphy 

of  poriaf  WjUdO  k  is  —064*  to  exploit  optical 
prnn— g  m  ikt  — i i  pkiingr— r  sugt  to  fttrtkar 

t—^pkoMfrapkic  ayMOMaai  kotgnd  onkipit  (ia- 
dodiof  iMy)  of  a  §wta|  frapraacy.  Sock  i  proodov 
«  oaafid  for  ■■kiptyit  aMonvky  or  for  nkttaaag  km- 
—  rttdok|  of  poii|i  aod  fnapt  potm.  Drptk 
of  ffaM  k  Ircrwidg  aod  a  coon  law  of  poor*  —key 
tkoo  io^aonaaky  r— aad  far^JOp^wokokaa  pkoto- 

Ckn^aaddoafli— an«a— —  km.  ** 
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